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Abstract: We investigate the problem of prolonging network life time considering an upper bound on the tolerable delay. 
According to the characteristics of the application in-hand, we put a constraint on the number of active nodes. In the 
proposed approach a mobile sink is taken advantage. The problem is formulated as a Linear Programming optimization 
using which both sink's mobility pattern and the optimum flow assignment are obtained. Using the proposed formulation, a 
(1-e) optimal solution (i.e., the network lifetime is at least (1-s) of the optimal network lifetime) is obtained. Numerical 
results indicates that by jointly considering the QoS requirements, routing, scheduling and sink mobility the network 
lifetime can be improved significantly. 
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I. Introduction 

Wireless sensor networks (WSNs) consisting of a large 
number of wireless sensor nodes which are spatially 
distributed in a designated region. Sensor nodes provide a 
distributed environment sensing for remote or inaccessible 
areas. In fact, each sensor node can gather information from 
its surrounding area. 

The collected data will be delivered to the base station, 
which is also known as a sink node, generally via multi-hop 
communication. The sensor nodes are typically expected to 
operate with batteries while battery replacement is an 
inconvenient and expensive task and usually is not possible. 
As a consequence, efficient energy utilization for prolonging 
the network lifetime has been the focus of much of the 
researches on the WSNs [1]. 

Sensors are typically deployed densely in the monitoring 
region. Due to spatial redundancy of sensor nodes and wide 
sensing range, same part of the monitoring field may be 
sensed by different sensors simultaneously. 

On one hand, in critical application such as battlefield 
surveillance, overlapped monitoring makes the collected 
data more reliable and because of that all sensors are 
required to be in active move. 

On the other hand, in some applications where network 
lifetime prolonging is more critical than providing accurate 
information, network lifetime can be prolonged by using 
only a subset of sensors at each point of time. For instance, 
in temperature sensing applications, sensing the temperature 



of a point could be adequate to its surrounding area. Thus, 
putting a subset of sensors into sleep mode to save energy is 
an efficient approach for network lifetime elongation in 
these kinds of applications. 

Developing a real-time applications over WSNs should 
consider not only resource constrains, but also providing a 
certain delay guarantee for data delivery. 

In contrast with the real-time applications, in delay- 
tolerant applications, employing data buffering methods is 
an efficient way for prolonging network lifetime. In other 
words, in applications which can tolerate a degree of delay 
in delivering the data packets to the sink node, each sensor 
node can postponed its data transmission for a while. 
Especially when a mobile data gathering is used, each sensor 
node can postponed its data transmission until the sink node 
is closed enough. 

In fact, WSNs are quite useful in many applications such 
as battlefield surveillance, medical diagnostics, 
environmental and habitat monitoring, each having specific 
requirements. In order to deliver the collected data with 
different requirements in highly energy constrained WSNs, 
application QoS provisioning becomes unavoidable. 
Actually, QoS support mechanisms brings the ability of 
giving different priorities to varied applications according to 
their requirements. 

QoS provisioning are utilized to manage the resource 
sharing according to the application-specific requirements. 
For the purpose of application differentiation, we focus on 
two QoS parameters listed as follow: 
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1- Accuracy, what percent of sensor nodes are needed 
to be in active mode at each point of time to meet the 
application accuracy level? 

2- Upper-bound on tolerable delay: How much data 
delivery delay is tolerable based on the application in-hand? 

The question is to find a routing paths meeting the 
minimum application-specific QoS requirements. In 
practice, longevity of the network lifetime is desired in many 
applications. 

Due to the limited battery power of the sensor nodes, it is 
extremely important that the routing be energy efficient, 
which aims at increasing the network lifetime. Routing 
protocols aiming to save energy in WSNs mainly focus on 
the sensor nodes, whereas a recent trend indicates a focus 
shift to the behavior of sink node, which can be employed to 
further improve the lifetime of WSNs. In fact, using a 
mobile sink or relay node can improve the network lifetime 
significantly [2-4]. 

The main reason for the improvement brought by the sink 
mobility approaches is the typical many-to-one traffic 
pattern in WSNs which is also known as convergecast. Such 
a flow pattern imposes a heavy forwarding load on the nodes 
close to the sink node. As a result, sensors which are closer 
to the sink node tend to run out of energy faster than the 
other ones. 

While no energy-conserving protocol alleviates such a 
load, moving the sink node can distribute the role of 
bottleneck nodes over time. Furthermore, mobile sinks may 
improve the network connectivity by enabling the retrieval 
from several isolated parts of the WSNs. 

In this paper, we exploit sink mobility to maximize the 
network lifetime for a WSN with a known application QoS 
requirements. For (1 - e) optimality, the infinite possible 
locations for sink mobility reduced to a finite set of 
locations. 

Our main contribution is that the proposed formulation 
fulfill the application QoS requirements, and in addition, it 
results a (1 - e) optimal solution. The problem is formulated 
based on linear programming optimization subjected to 
maximize the network lifetime while meeting application 
QoS requirements are guaranteed. 

The rest of this paper is organized as follows: In Section 
II, we have given a brief overview of the state-of-the-art on 
related lifetime maximization problems. In Section III, we 
gave a problem definition and statement in detail. Proposed 
method is presented in section IV. Performance evaluation 
and simulation results are shown in section V. In section VI, 
the paper ends with a concluding epilogue along with a hint 
on future works conceivable in this area. 



II. Related Works 

In this section, we present a comprehensive review on the 
related lifetime maximization problems that have been 
published in the literature. 

Optimal sink node placement has been an active area of 
research for network lifetime maximization in many recent 
studies (see, e.g, [1-3]). The focus was mainly devoted to 
determination of the optimal sink location and sensor-to-sink 
information flow routes. However, network lifetime can be 
further improved by using a mobile sink [2-6]. 

In [6], Wanget et al tried to maximize the network 
lifetime by considering a pre-determined routing paths. In 
fact, they proposed a new formulation for determining the 
sink sojourn time in order to maximize the network lifetime. 
This differs from our proposed method in that we jointly 
consider routing, QoS and sink mobility. 

However, in [8] Papadimitriou and Georgiadis extend the 
formulation used in [6] by jointly considering sink mobility 
and routing to maximize the network lifetime. The sink node 
is however confined to a limiter set of location (i.e., the sink 
node can just move between a predefined set of anchor 
points). 

In [7] authors investigated on the joint sink mobility and 
routing to improve the network lifetime. Limited set of 
locations are used as the sink anchor points. They proved the 
superiority of mobile sink node over static sink node in the 
case that the sinks are constrained to where the nodes are 
(i.e., the set of sensor locations selected as possible anchor 
points for the sink mobility). They proposed a linear 
programming formulation in which the sink sojourn time 
and corresponding routing are determined in order to 
maximize the network lifetime. 

However, they did not considered the transmission 
distance in the transmission energy consumption model 
whereas the recent studies showed that tuning the 
transmission power according the transmission distance can 
improve the network lifetime [10]. 

Besides, their formulation is not applicable for 
unconstrained sink mobility. An extension of [7] is presented 
recently in [9] where an approximation method used to 
develop the constrained sink mobility into unconstrained 
sink mobility. However, they did not assumed QoS 
requirements in their formulation. 

Distributed sink mobility and routing approach proposed 
in [11] in which the sink node can move between a limited 
and predefined set of locations. They proposed a linear 
programming formulation for joint sink mobility and routing 
in order to maximize the network lifetime. Lagrangian [19] 
and subgradient [20] method used to solve the problem in 
distributed manner. However, they did not consider the 
energy consumed in a sensor node for data reception. Their 
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proposed method consists of two different phases. 
Initialization phase requires co-operation between sensor 
nodes to achieve a near optimal solution. Besides, in routing 
phase, the solution found in the initialization phase is used to 
forward the data packets toward the sink node. However, 
low convergence speed of lagrangian and subgradient 
method cause energy depletion of some sensor in the 
initialization phase. 

Same as [11], in another contribution [12] authors 
investigated on the distributed algorithm for mobile data 
gathering in wireless sensor networks. They employed 
anchor based mobile data gathering to maximize the network 
utility under the constrained and guaranteed network 
lifetime. 

In fact, the only parameter assumed as QoS is the 
network lifetime which is guaranteed to achieved while the 
network utility is maximized. They used Lagrangian 
heuristic along with subgradient method to solve the 
problem in distributed manner. However, they did not 
assumed energy consumed in a sensor node for data 
reception. Besides, their proposed method has low 
convergence speed. 

There were some efforts on lifetime maximization 
problem by using mobile sink in delay tolerant applications. 
[13-15]. Yun and Xia[13] exploited the sink mobility in 
WSNs where undelaying applications tolerate delayed 
information delivery to the sink node. 

They proposed a new formulation in which each sensor 
node sends data to the sink when the mobile sink is within a 
given distance of the sensor. In fact, at each point of time, 
only sensors which their distance to the sink node is less or 
equal than specific value should send their data to the sink 
node. Obviously, using this approach cannot guarantee the 
end to end delay in data delivery. In contrast with [13], in 
our proposed approach, end to end delay is guaranteed to be 
less or equal than specific value. 

In [14] authors studied on joint optimization problem 
consisting of trajectory finding and sojourn time scheduling 
with an aim to maximize the network throughput. They 
assumed a set of potential sojourn locations for the mobile 
sink node and employed a heuristic method to find out the 
sink sojourn time and routing at each location. 

III. System Model 

We model a WSN as digraph G = (V,E) where \V\=n 
indicates the sensor nodes and 3(t,j) EE if and only if 
node j is located within the transmission range of node i. 

Let S represent the sink node which harvest data from the 
WSN. All sensor nodes are stationary while the sink node 
can change its location from time to time with a negligible 



traveling time between two locations. Each sensor node t 
which is located at {xi.yi) has E Joules of energy as initial 
energy and generates data with the rate of r. 

The generated data by each sensor node should be 
transmitted to the sink node S via a single hop or a multi-hop 
communication. Table I lists the notation used in this paper. 

We assumed that each sensor node can control its 
transmission range. 

Suppose node i and j while d shows the Euclidean 
distance between node i and j. Let E-j denotes the energy 
consumed for sending a bit of data from node i to node j. 
Then the transmission energy cost is modeled by 

E!j=e t + I3.d?j (1) 



Where e t and /? are constant coefficients and a is path 
loss parameter which takes value between 2 and 4. 

Let E r denotes the energy consumed for receiving a bit of 
data. Since the sender node controls its transmission power 
level, the receiver node consumes a constant amount of 
energy for receiving a bit of data. 

We assumed that the sink node can change its location 
from time to time with negligible traveling time. Since a 
mobile sink node is employed, in delay-tolerance 
applications each sensor node can postpone its data 
transmission until the sink node stops at the most favorable 
location for extending the network lifetime. 

Each sensor senses its surrounding area periodically and 
generates data packets with a certain rate. Each sensor has 
limited amount of memory which can be employed for 
buffering data packets. In other words, each sensor is able to 
gather data form its surrounding environment and buffer 
data by using its limited memory and send them to the sink 
node in a certain period of times. 

Based on the application, each node can postpone its data 
transmission for a certain amount of time. Let D indicates 
the maximum tolerable delay in seconds. 

By using a FIFO queuing policy for selecting data packets 
from the memory buffer, for sensor node i which generates 
data with the rate of r bits per second and has M bits of 
memory for data buffering, maximum data delay would be — 
second. 

Note that we ignore transmission and propagation delay. 
Besides, we assumed that each sensor node use its buffer 
memory to store its own data packet. 

As a consequence, by controlling the buffer queue length, 
one can insure end to end data delivery delay 
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Table I. notations 







s 


Sinknode S ° rS 


£ 


Set of wireless links 


a 


Path loss parameter 


E h 


Energy consumed for sending a bit of data 
from sensor ;' to sensor j. 




Energy consumed for receiving a bit of 
data. 




Sensor initial energy. 


r 


Data generation rate 


A 


Center of SED 


Ra 


Radius of SED 


D min 


Minimum distance from sensor node i to 
the sink node 


D rnax 


Maximum distance from sensor node i to 
the sink node 


D io A 


Euclidean distance from sensor node ;' to 
the center of SED 


E_min\ s 


Minimum Energy required for sending a bit 
of data from sensor node ;' to the sink node 


E_maxf s 


Maximum Energy required for sending a 
bit of data from sensor node i to the sink 
node 




Transmission energy cost for sending a bit 
which is located at subarea j 




sensor node i to cover SED completely. 






A m 


Subarea m 


Vm 


Point located at A m 


T 


Network lifetime 


D 


Endurable delay 


P 


Accuracy parameter 


F ij 


Amount of data flow send by node / to 
node j during k-th epoch. 




Amount of data which node stores at its 
own data buffer during k-th epoch 


Boutput i 


Amount of data which node ;' sends from its 
own data buffer during k-th epoch 


Blength* 


Queue buffer length of sensor i during k-th 
epoch 


Blength i 


Queue buffer length of sensor i at the 
system initialization 



We now argue on how to move the sink node with (1 - e) 
optimal network lifetime guarantee. Denote A the movement 
region for the sink node, which can be narrowed down to the 
smallest enclosing disk (SED) for all sensor nodes, i.e., 
smallest circular area that covers all the nodes in the network 
[15]. Note that for a given WSN, the SED can be found in 
polynomial time [16]. Denote A and R A as the origin and 
radius of SED A respectively. Let D ™ in and D ™ ax the 
minimum and maximum Euclidean distance from sensor i to 
sink node S. Since the movement region for the sink node S 
is within SED A, we have: 



DT = (2) 
D% ax = D i0A + R A (3) 

Corresponding to D ™ n and D™ ax , denote Ejnin^ and 
E_maxf s the minimum and maximum transmission energy 
cost between sensor node t to the sink node S. Then by 
(1) ,(2) and (3), we have: 

Ejnin\ s = e t (4) 
Emaxls = e t + KD^ ax ) a = e t + B{D Wa + R A ) a (5) 

Obviously, the Euclidean distance between sender and 
receiver affects the transmission energy cost. Thus, to 
achieve (1 - e) optimal network lifetime, we divide SED 
into finite number of subareas while the Euclidean distance 
from each subarea to sensor i meeting some tight bounds. 
For each sensor £ V , each subarea j is associate with 
transmission energy cost Ef s \j] which shows the 
transmission energy cost for sending a bit of data to the sink 
node 5 which is located at the subarea j. 

In fact, discrete energy cost is used based on geometric 
sequence [18]. Specifically, by drawing a sequence of circles 
centered at node t, each with increasing distance 
d[2], d[Hi] corresponding to transmission energy 
cost Ef s [l],Ef s [2], ... ,Ef s [Hi] which are defined based on 
the following geometric sequence: 

Ef s [h] = E_min\ s {l + e) a = e t (l + e) a 1 < h < H i (6) 

In fact, E\ S \\\,E\ S \2\, ... ,E\ S \H?[ is geometric sequence 
with the factor of (1+e). 

For each sensor i, the number of required circles to cover 
A completely (i.e., H_i) can be determined by having the last 
circle in the sequence which : 

£&[//;] > Emaxfs (7) 
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Thus Hi can be determined as follows: 

_ UnCEjnaxfs/Ejninjs)] 
Hi ~ Tn(l + e) | 

_ ln(l+§(D i0A + R A ) a ) ] _ /Mix 
In(l+e) ~ \\s\) 

= 0y (8) 

Where the third one hold by ln(l + e) « e for small e and 
1 + ^7 (Pw A + Ra)" is a constant. For each sensor node i, 
drawing the circles according to the above procedure results 
Hi non-overlapping rings each has radios of < j < 

H t ) . By locating the base station between (j-1 J-th circle and 
j-th circle, the transmission energy cost for sending a bit of 
data from sensor node i to sink node S would be: 

Ef s [h - 1] < Efs < Ef s [h] (9) 



While we assume that E- s [0] = E_min\ s = e t . Because of 
f'f , = 1 + £ based on (6), we can claim that these two 

E is [h-1] 

bound for Ef s are very tigh. 

By performing the above procedures for all sensor nodes, 
the intersecting circles will divide disk A into a finite 
number of subareas with irregular shapes. Note that for each 
sensor node i, any subarea A m is located within circle 
centered at node i. Denote the index of this circle as 
h t (A m ) . So, for a sink node S which is located at any point 
within A m , transmission energy cost would be: 

Efs < Ef s [hi(A m )] (10) 



So each subarea A m assumed as point p m where the 
transmission energy cost for sending a bit of data from 
sensor to the sink node S located at p m assume as follows: 

El; = Et s [hi(A m )] (11) 



For unconstrained sink mobility model we can assume the 
SED as a set of finite points. 

The question is to find a sink sojourn time at each point 
p m and corresponding data routing solution which maximize 
the network lifetime. 



IV. Problem Statement and Formulation 

We assume a set of sensors which distributed randomly in 
a monitoring region. We assume that the location of sensor 
nodes and sink nodes are known at the system startup. If cat 
we can used GPS based or GPS less localization method to 
achieve the sensors location. 

Sensor nodes gather data from monitoring region while 
the sink node changes its location from time to time in order 
to gather data from sensor nodes. 

Based on the predefined QoS requirements, at least p 
percent of all sensors should generate data with the rate of r 
bit per second at each point of time while each sensor node 
can postpone its data transmission for D seconds 

Based on the sink mobility model, the sink node can 
move around the SED while the SED is divided into a finite 
number of subareas. 

Since the transmission energy cost for each subarea is 
constant for each sensor node, we assume each subarea as 
point. The sink node can stop at each point for a certain 
amount of time. 

Although the network lifetime can be defined in many 
ways, we adopt the time until the first sensor node exhausts 
its energy as network lifetime, which is a widely used. 

We denote the network lifetime by T and use t k to 
indicate timespan the for the k-th epoch. By changing the 
sink node location, a new epoch began. 

We denote the total information flow from sensor node ;' 
to sensor node j during the k-th epoch by F-j. Based on the 
problem definition, for a given value of p, at least p% of all 
sensors should generate data by the rate of r at each point of 
time. Since we assumed that based on the application in- 
hand, not all sensors are needed to be in active mode, q' ik 
used to indicate total information flow which sensor does 
not send generate during the fe-the epoch. 

However, for delay-tolerance application, the memory 
buffer can be employed for storing data packets. 

Let B input ^ indicated the total information flow which 
sensor stores at its own buffer memory during the k-th epoch 
and B output k as total information flow which sensor sends 
from its own buffer during the k-th epoch. 

Without loss of generality, we assume that the sink node 
dwells at the k-th subarea during the k-th epoch. Note that 
this assumption does not change the problem. It only 
requires that the number of subareas coincide with that of 
epoch. 

Based on the flow conservation rule, the outgoing flow 
exceed the incoming flow by the amount of r X t k . 

So the flow conservation rule for each sensor node i at 
each epoch can be formulated as follows: 
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V;6VUS 



Vjev 



Vi.fc 



Ufa 
(12) 



vi (l8) 



vfc v/ev vk vjev 



Where E r is the energy consumed for receiving a bit of 

Obviously, each sensor node can sends its data to another data and Efj is the energy consumed for sending a bit of data 

sensor node, sink node or it can store in its own data buffer from sensor node f to node j defined by (1). Besides, 

memory. energy consumed for sending a bit of data by sensor node i 

Let B length k . be the buffer queue length of node / during to the sink node during the k-th epoch defined by (11). 
the k-th epoch. Since the buffer is empty at the system As a result, we can formulate the problem as linear 

initialization we have programming formulation as follows: 



a length i 



~- 



(13) 



Max T = 



(19) 



During each epoch queue length is the sum of queue 
length in the previous epoch and data generated but not 
transmitted in that epoch. As a result, the queue length at 
each epoch can be formulated as follows: 

u*-B 0Utput * Vi.fc (14) 



(12) - (18) 



t >iength i : 



° length t 



In which during the k-th epoch which last for t k seconds, 
sensor ;' generates data with the rate of r bits per seconds and 
transmit q ik of the generated data. Sensor node i buffers Bf 
bits of data during the k-th epoch .Since we assume that each 
application has a specific data-tolerance D, the queue length 
at each sensor should not exceed D. i.e.: 



o length i 



Vi.fc 



(15) 



Besides, we assume that based on the predefined QoS 
requirements, at each point of time p percent of all sensors 
should generates data. For k-th epoch which last for 
t fc seconds, at least p X \V\ of all sensors should generate 
data with the rate of r bits per seconds, i.e.: 

i(l-p)\V\rt k Vfc (16) 



Besides, at the end of the last epoch, all sensors should 
empty their buffer and send all data to the sink node. In other 
words, for the WSN with T seconds of network lifetime and 
| V | sensor nodes, at least p% of all sensors should generate 
and transmit data with the rate of r for T seconds. 

II 



Z(l-p)\V\rT 



(17) 



However, each sensor has limited initial energy denoted 
by E . Total energy consumed for data transmission and 
reception should not exceed the initial energy. In other 
words: 



The outputs of the above formulations are the network 
lifetime, sink sojourn time at each subarea, flow assignment 
for each sensor node at each time epoch, amount of data sent 
by each node during each time epoch, amount of data which 
does not send by each sensor node during each time epoch, 
queue length for each sensor at each time epoch. 

V. Correctness Proof and Complexity Analysis 

In this section, we give a formal proof that the solution 
obtained by Linear programming formulation is (1 - e) 
optimal and analyze its complexity. 

Denote {pi, p 2 , ... , p n } as the optimal anchor points for the 
sink mobility, {t£, t 2 , ... , t*} and {R{,R 2 , ... ,R* n } as their 
corresponding sojourn times and routing respectively which 
are all unknown. Beside, let {p 1 ,p 2 , ■■■ ,Vn\ indicates the 
optimal anchor points (selected subareas) resulted by solving 
the linear programming formulation and let {t ± , t 2 , ... , t n } 
and {R 1 ,R 2 , ... ,R n }denoted the corresponding sink sojourn 
time and routing respectively. We also denote the T* and T 
as optimum network lifetime and network lifetime resulted 
by solving the proposed LP. 

Theorem. Based on the definition of T* and T , we have 
T > (1 - £)7* 

Proof. At each time epoch k, instead of using the optimal 
routing R k we use R* k for the solution resulted by LP which 
is clearly suboptimal. That is, denoting t as the network 
lifetime for LP under R* k for (1 - e)T, we have T > f . Then 
it is only required to show that t > (1 — s)T* 

Since we use each R* k for at least (1 — £)t- instead of tj 
seconds, for each £ R k and F*fj £ R k we have 
(1 - £)F% <F t k j< F*1j (20) 
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Then to prove f > (1 - e)T* , we only need to show the 
total consumed energy at each sensor node does not exceed 
its initial energy, which is: 



The first inequality holds by (6) . Thus the network 
lifetime for the solution of the LP under R* is at least 
(1 — a)T*. As a consequence, we have (1 — e)T* <f< 
T <T* which completes the proof.H 

The number of subareas has great impact on the problem 
complexity and the number of LPs which needed to be solve. 
Thus to calculate the problem complexity, we are going to 
approximate the number of subareas. 

The boundaries of each subareas is either an arc of the 
circle centered at some node t or an arc of the disk A. Since 
for each sensor node there are H t — 1 (defined by(8) ) 
circles and there is a circle for A there are totally N = 1 + 
T.viEv(.Hi — 1) cicles. Based on [20] the total number of 
subareas M can be upper bounded by: 

M<N 2 -N + 1 (21) 

Based on (8) and (21) we have 

M = 0(JV 2 ) = 0((JV/ £ ) 2 ) (22) 

It can be easily seen that the number of constraints and 
variables used in the proposed LP is about 0(|7| 2 ) and 
(|K| 3 ) . Hence, the problem can be solved in polynomial 
time [21]. 



VI. Performance study 

In this section, we evaluate the proposed method using an 
extensive numerical study. The parameter values used 
during the numerical study are listed in table I. To achieve 
0.95 optimal solution we set e= 0.05 . We investigate three 
metrics, namely network lifetime, sink sojourn time 
distribution and sensors buffer queue length. In order to 
study our approach for different network topologies we take 
into account three types of network topologies including 
samples of linear, grid and arbitrary topologies as follows. 
The symbols definitions are accessible in table I. 

The proposed LP is compared with the ones proposed in 
[10] and [6]. In fact, in [6] the authors proposed a sink 
placement approach based on a mixed integer linear 
programming formulation. Since the integer programming 
formulation is NP-hard, they proposed an efficient heuristic 
to find out the near optimal sink location. Meanwhile, in 
[10] authors investigated on joint sink mobility and routing 



while the set of anchor points for the sink node are 
constrained to the sensors locations. However, neither [6] 
nor [10] assumed any sort of QoS provisioning. However, in 
our approach we consider an upper bound on delay as the 
QoS requirement. 



Table II. 



Parameter 


Value 


r 


50bps 


a 


3 


e t 


50 nj/bit 


E r 


150 nj/bit 


£ 


0.05 




40 m 


E 


5 Jules 



In the line network shown in Fig.l, the distance between 
each two neighborhoods are selected within [20,40]. 



O- -O -O O O O O- O -O 

Fig. 1 A line network topology 



As shown in Table III, using unconstrained mobility (in 
our approach) can always achieve a higher network lifetime 
than using a static sink node. 

Specially by increasing the size of network, the 
improvement in network lifetime is more significant. 
Comparing with the constrained sink mobility ( in which the 
candidate locations for the sink node include only sensor 
locations) the proposed approach ( in which candidate 
locations for the sink are some discretized subareas) 
improves the network lifetime. 

Besides, it can be seen that if the application can tolerate 
200 seconds (i.e., D=200) of delay in data delivery the 
network lifetime improves about 1.8 in comparison with 
unconstrained mobile sink with D=0. 

Furthermore, by putting only 20% ( i.e., p=80%) of 
sensors into sleep mode the network lifetime improves about 
2.09 in comparison with the case that p = 100% . 

An interesting point about the line network topology 
which is also confirmed in our study and previously 
observed in [10] is that the improvement network lifetime is 
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mostly resulted by the substitution effect [10]. That is, the As a result, sensors do not gather data from their 

sink node inherits data forwarding load from the neighbor surrounding area with a same rate in these kinds of 

nodes. applications. Fig.2 shows the normalized rate of not 

As we will show in the next section, for other topologies, generated data while D=0 and p=80% for N=49. The sink 

the lifetime improvement is caused mostly by load balancing node tends to dwells mostly in the center of the target field, 

effect rather than substitution effect of mobile sink. Therefore, sensor nodes which are far from the sink node 

For grid network with J\v~\ X J\v~\ lattices, the average S ather data from their surrounding are with a lower rate than 

network lifetime achieved for 100 random generated sensors which are located in the center of the field, 

topologies are shown in Tables III We a ^ so consider network with N=49,D=200 and p=80% . 

As followed from the table, for large scale networks, Fi S- 2 shows the sensors buffer queue len S th of such a 

using the static sink node cause unbalance energy network. Considering Fig.2 and Fig.3 shows that sensors 

consumption resulting a short network life time. However, which far from me sink node (i- e - mostl y located around 

as the network size increases, the improvement on the the tar S et field ) tends to generate data with a lower rate 

network lifetime resulted from unconstrained sink mobility while me nearb y ones have more data in their buffer queue, 

increases. We a l so P er f° rm the proposed method on arbitrary 

Using a mobile sink node in grid network improves the network topology (i.e., sensor nodes are deployed based on 

network lifetime not only by the load balancing effect but the uniform distribution on the target field), 

also by the substitution effect. Because of that, the Fi S- 4 shows the avera ge network lifetime over 100 

unconstrained sink mobility improves the network lifetime random generated topologies for different scenarios while 

about 2.3 times while D=0 and = 100% . network size varies from 50 tol20 and D=0. Clearly, the 

However, in delay tolerant applications when 200 seconds network lifetime improves by decreasing the level of 
of delay is tolerable using the mobile sink improves the ^curacy- 
network lifetime about 2.6. Besides, using at least 80% of In critical plication with high accuracy requirement one 
sensors at each point of time improves the network lifetime can use = 100% . However, using p = 90% and p = 80% 
by 3.1. Furthermore, when 200 seconds of delay in data improve the network lifetime by about 16% and 28% 
delivery is tolerable and at least 80% of all sensors are active respectively. 

at each point of time the network lifetime improves by 3.78. Fi S- 5 illustrates the comparison of network lifetime 

In some applications the required accuracy level can be achieved by different scenarios. For real-time applications 

fulfilled by using only a part of sensor nodes instead of using we can use D = 0. However, when the D = 200 seconds is 

all sensors. Thus, only subsets of sensors are selected as tolerable by application the network lifetime improve by 

active sensors at each point of time. °- 25 - Moreover, using D = 400 and D = 600 improves the 

network lifetime by about 0.35 and 0.55 respectively. 

Table III Network lifetime resulted by different methods in linear topology 

Network Lifetime 

Proposed method ( Unconstrained mobile sink) 

W\ 

Static sink Constrained 

placement mobile sink (D=0,p=100%) (D=200,p=100%) (D=0,p=80%) (D=200,p=80%) 
[6] [10] 

20 345.76 386.20 458.12 624.72 723.76 918.82 

40 726.096 811.02 958.4467 1309.151 1437.67 1103.666 
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Fig. 3 Average Buffer queue length for a network with N=49, 
D=0 and p = 80% 




Fig.4 Impact of network size and level of accuracy on the 
network lifetime 



VII. Conclusion 

In this paper, we present a new framework to achieve 
the maximum network lifetime for different kinds of 
application with particular QoS requirements. We exploit 
joint sink mobility and routing to maximize the network 
lifetime while QoS requirements are considered as well. 
For the purpose of QoS provisioning, we consider two 
metrics: 

a) Accuracy level: what percent of sensor nodes are 
required to be in active mode at each point of time based 
on the application in-hand. 

b) Handling the delay: How much data delivery to the 
sink node is tolerable? In delay tolerance application each 
sensor node can postpone its data transmission until the 
sink node is closed enough. However, in uncritical 
application with low accuracy level requirement, only a 
subset of sensors could be used to improve the network 
lifetime. 

For sink mobility model, we used an unconstrained 
sink mobility to achieve (1 — s) optimal solution. In fact, 
for (1 — £) optimality, the target region can be 
discretized into a set of subareas, each having a specific 
transmission cost. Problem formulated based on the 
linear programming formulation with the respect of 
predefined QoS requirements. 

Besides, the correctness and complexity analysis 
presented which proves the (1 — £) optimality of the 
solution resulted by solving the proposed LP. More 
important, we prove that the proposed solution can be 
solved in polynomial time. 
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